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Abstract: Hard exclusive leptoproduction of real photons from an unpolarized proton 
target is studied in an effort to elucidate generalized parton distributions. The data ac- 
cumulated during the years 1996-2005 with the HERMES spectrometer are analyzed to 
yield asymmetries with respect to the combined dependence of the cross section on beam 
helicity and charge, thereby revealing previously unseparated contributions from deeply 
virtual Compton scattering and its interference with the Bethe-Heitler process. The inte- 
grated luminosity is sufficient to show correlated dependences on two kinematic variables, 
and provides the most precise determination of the dependence on only the beam charge. 
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1. Introduction 



Parton Distribution Functions (PDFs) describe the longitudinal-momentum structure of 
the nucleon in the interpretation of inclusive and semi-inclusive Deep-Inelastic Scattering 
(DIS). Analogously, Generalized Parton Distributions (GPDs) ||l],|2|, ^ describe the multidi- 
mensional structure of the nucleon in the interpretation of hard exclusive leptoproduction, 
most simply when the target is left intact. PDFs and elastic nucleon Form Factors (FFs) 
are embodied in GPDs as their limiting cases and moments, respectively |^. While FFs 
and PDFs represent one-dimensional distributions, GPDs provide correlated information 
on transverse spatial and longitudinal momentum distributions of partons [Q, |5|, ^. In 
addition, the total angular momentum carried by partons in the nucleon can be calculated 
from GPDs @. 

GPDs depend on the four kinematic variables x, ^, Q"^, and the squared four- momentum 
transfer t to the target. In a frame where the nucleon moves with 'infinite' momentum, x 
and 2^ are the average and difference of the longitudinal momentum fractions of the parton 
in the initial and final state, as illustrated in Fig. 1(a). In hard exclusive leptoproduction, 



X has no direct relationship with the experimental kinematic observable xb = Q^/{2P ■ q). 
Here, P is the four momentum of the target nucleon, q is the difference between the four 
momenta of the incident and scattered lepton, and = —q^- The skewness parameter ^ 
is related to xb, as ^ xb/{2 — xb) in leading order Quantum Electrodynamics (QED) 
and in the kinematic limit of large with xb and t fixed (generalized Bjorken limit). 
In addition, like PDFs, GPDs are subject to Quantum Chromodynamics (QCD) evolu- 
tion with Q^, which has been calculated perturbatively to leading order |2|, |3|, |^ and 
next-to- leading order |8|, |9|, ^ in the strong coupling constant as- 





(a) 



(b) 



Figure 1: Leading order diagrams for deeply virtual Compton scattering (a) and Bethe-Heitler 
(b) processes. 

Deeply Virtual Gompton Scattering (DVGS), the hard exclusive leptoproduction of real 
photons, e.g., e^p — > e^p7, has the simplest theoretical interpretation in terms of GPDs 
among the presently experimentally feasible hard exclusive reactions. DVGS amplitudes 
can be measured through the interference between the DVGS and Bethe-Heitler (BH) 
processes, in which the photon is radiated from a parton in the former and from the 
lepton in the latter (see Fig. |l|). These processes have an identical final state. Hence their 



-3- 



amplitudes tdvcs s-nd tbh add coherently, resulting in an interference term ' I ' in the cross 
section for exclusive leptoproduction of real photons. For an unpolarized proton target, 



the cross section can be written as [|ll|, [12 

dfj xb 
dxB dQ^d\t\d(j) ~ 32 (27r)4 Q4 ^/f^ 



ItbhP + koVCSp + ''"DVCS ^BH + "^DVCS ''"BH 



a.i) 



where e represents the elementary charge and e = 2xBMp/Q, with Mp the proton mass. 
The azimuthal angle (p is defined as the angle between the lepton scattering plane and the 
photon production plane spanned by the trajectories of the virtual and real photons. 

The three contributions entering the photon production cross section can be expanded 



in Fourier series in (/). For an unpolarized proton target, they can be written as |12 



\^ |2 

rBHi 



PDVCSl 



KbH BH / ,N I ■,, 1 



^c^^cos{n(l))\ , with ^bh = ^— 7— (1-2) 



^ IEc^v^'cosM) + A^r^^sin^l , (1.3) 
cjj cos(n0) + A sin(n</)) > , with Ki = . (1.4) 



^l(</')^2(0) " ^1 " J ^B2/t 

Here, y is the fraction of the incident lepton energy carried by the virtual photon in the 
target rest frame, and A and represent respectively the beam helicity and beam charge in 
units of the elementary charge. The Fourier coefficients and lepton propagators Vi{(j)), 
V2{4>) of the BH term can be calculated within the framework of QED from the kinematic 
variables and the Dirac and Pauli form factors -Fi and F2 of the nucleon. 

The Fourier coefficients of the interference term (Eq. |1.4| ) are of greatest interest since 
they ultimately depend on a linear combination of GPDs, while the coefficients of the 
squared DVCS term (Eq. |1.3| ) are bilinear in GPDs. The coefficients 

s\ = 8kXy{2-y)lmM^^\ (1.5) 
c\ = 8k{2 -2y + y^) ReM^^'^, (1.6) 

are respectively proportional to the imaginary and real parts of M^'^, the leading-twist 
(twist-2) photon-helicity-conserving amplitude of the DVCS process. Here, M^'^' denotes 
the helicity amplitude with virtual (real) photon helicity /x (/i'), following the notation of 



Ref. 1 11]. The kinematic factor k oc ^/—i/Q originates from the BH propagators. Note 



that the sign in Eq. |1.6| differs from that in Ref. |1^ due to the different definition of the 

IT2I 



azimuthal angle: (j) = tt — (pffp^- The amplitude M^'^ is given by a linear combination of 



the Compton Form Factors (CFFs) 7i, 7i and £: 
Mi'i = F^it) m^, t, Q2) + (F, (t) + F2it)) n{^, t, Q^) - ^F^it) £{i, t, Q^). (1.7) 
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The CFFs are convolutions of the respective twist-2 GPDs H, H or with perturbatively 
calculable hard-scattering amplitudes. These amplitudes have been calculated to next-to- 
leading order in as 14, 15 1. The contributions from the CFFs Ti and E to the amplitude 



M^'^ are kinematically suppressed compared to that from the CFF TL at small values of 
xb and t, respectively. 

In addition to s\ and c^, the only other Fourier coefficients related to quark-helicity- 
conserving twist-2 GPDs are Cg and Cq^*""^. The coefficient Cq is also related to M^'^. 
Considering only the dominant CFF Ti^ Cq is directly proportional to (\ via the factor k 
defined above: 



Co oc 



-fcc^. (1.^ 



The coefficients C2, s^^*""^ and cj^^^^ are related to twist-3 GPDs. The coefficient s\ 
(C2) is proportional to the imaginary (real) part of the helicity non-conserving amplitudes 
M^'^^ , corresponding to the virtual photon being longitudinal. Conservation of angular 
momentum is ensured by either the exchange of an additional gluon (genuine or dynamic 
twist-3) or by the fact that quarks can carry non-zero orbital angular momentum along the 
collision axis (kinematically suppressed by the same order in 1/Q), which is possible due 
to the transverse momentum involved. The part of the twist-3 GPD associated with the 
latter picture can be related to the twist-2 quark GPDs using the Wandzura-Wilczek (WW) 
approximation [^] and thus is also known as the WW part of the twist-3 contribution |1T7|| . 

The Fourier coefficient Cg is proportional to the real part of the amplitudes M^~^ and 
M~^'^, which do not conserve photon helicity, i.e., both photons are transverse and they 
have opposite helicity. The induced two units of angular momentum can be accommodated 
by gluon helicity- flip. Gluon helicity-flip GPDs do not mix with quark GPDs via 
evolution and thus probe the gluonic properties of the nucleon |10|. They appear at leading 
twist, but are suppressed by a factor Os/vr. In addition, as in the case of the coefficients 
discussed above that are kinematically suppressed by l/Q, it is possible that the two 
participating quarks complete the conservation of angular momentum if they carry orbital 
angular momentum. As they have to account for two units of angular momentum instead 
of one as above, this process appears at twist-4. The associated twist-4 GPD was found to 
be calculable in terms of twist-2 quark GPDs using the WW approximation [l^. Similarly, 
the Fourier coefficient c^^*-"^ arises from the twist-2 gluon helicity-flip GPDs with possible 
contributions from twist-4 quark GPDs. 



2. Asymmetries 



Previous measurements with a longitudinally (L) polarized positron {electron} beam by 
HERMES H] {CLAS H, ||, H} on an unpolarized (U) proton target provided access to 
a combination of s\ and sj^^^^ via the single-charge beam-helicity asymmetry, also denoted 
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as the Beam-Spin Asymmetry (BSA): 



-4lu(</', e^) = 



da 



da^ 



da~' + da'" 



Ki 


" 2 

J]] sj, sin(n0) 

.n=l 




1 


2 3 

-f^BH E c^^ cos(n(/)) - e^Ki ^ cj, cos(n(/)) 

n=0 n=0 


n=0 



(2.1) 



Here, a^ (o"^) denotes the cross section for a beam with positive (negative) hehcity. Pre- 
dominant sin (f) dependences with opposite sign have been observed at the two experiments, 
indicating the dominance of the interference term involving • s\. However, quantitative 
access to s\ is comphcated by the presence of s^^^^, which is a higher twist-contribution 
but possibly significant, and by the presence of (\ and Cg, i.e., the other Fourier coefficients 



of interest appearing at leading twist (see Eqs. |L6| and 

This entanglement can be avoided by defining the charge- difference beam-helicity 



asymmetry [23 



{da^ 



da^ 



{da- 



da 



{da+-* + da+^ 



+ {da — * + da-^) 

2 



n=l 



Vi(4>)'P2(4>) Z^"'" 

n=0 



E 



■DTT . 

c„ cosint? 



+ — Vc° 

^ n=0 



(2.2) 



DVCS 



cos nfi 



where the additional +(— ) superscript on the cross-sections denotes the charge of the lepton 
beam. This asymmetry has the important advantages that the siiKp dependence in the 
numerator stems solely from the interference term, as the (higher- twist) sin0 dependence 
of the squared DVCS term cancels, and the denominator no longer contains the leading 
terms c\ and Cg. Therefore it gives direct access to linear combinations of GPDs, while 
another charge- averaged asymmetry related to the squared DVCS term provides access to 
bilinear combinations of GPDs: 

(dfT+^ - (i(T+^) + {da — ' - da"-) 



{da+^ + da+^) + {da — ' + da- 



1 ^DVCS 



sm( 



(2.3) 



BH 



cos n<? 



+ 



n=0 



E 

n=0 



The previously extracted |24, 25 1 Beam-Charge Asymmetry (BCA) 

3 

[ cosine 



Ac{<t)) 



da^^ — da 
da^ + da- 



n=0 



2 2 

E ^n"" cos(n,^) + ^°^(™'^ 

n=0 ^ n=0 



(2.4) 
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provides access to the real part of the DVCS ampUtude via c]^. 
3. Event selection 



The data were cohected during the years 1996-2005 with the HERMES spectrometer |26| 
using the longitudinahy polarized 27.6 GeV electron and positron beams provided by the 
HERA accelerator facility at DESY. The hydrogen gas target was either unpolarized, lon- 
gitudinally or transversely nuclear-polarized. However, the time averaged polarization of 
the polarized targets was negligible, while the rapid (60 — 180 s) reversal of the polarization 
direction minimized polarization bias due to detector effects. The polarization direction 
of the beam was reversed about every two months. The integrated luminosity for the 
electron (positron) data sample corresponds to about 106pb^^(292pb~^) with an average 
magnitude of the beam polarization of 30.0% (49.4%). The latter has a mean fractional 
systematic uncertainty of 2.8 %. 

A brief description of the event selection is given here. More details can be found in 
Refs. 1 24, |2^. Events are selected with exactly one photon producing an energy deposition 



larger than 5 GeV (1 MeV) in the calorimeter (preshower detector) and one charged track, 
identified as the scattered lepton, in the kinematic range 1 GeV^ < Q2 < 10 GeV^, 0.03 < 
xb < 0.35, W > 3 GeV and < 22 GeV. Here, W denotes the invariant mass of the initial 
photon-nucleon system and u denotes the virtual-photon energy in the target rest frame. 
The angle between the real and the virtual photon is constrained to be between 5 and 
45 mrad. The recoiling proton is not detected. An 'exclusive' sample of events is selected 
by the requirement that the squared missing mass of the reaction p — > 7 X 
corresponds to the squared proton mass. The resolution in Mj- is limited by the energy 
resolution of the real photon in the calorimeter. Correspondingly, the exclusive region 
is defined as —(1.5 GeV)^ < M'^ < (1.7 GeV)^, as determined from signal-to-background 
studies using a Monte Carlo (MC) simulation. For elastic events (leaving the proton intact), 
the kinematic relationship between the energy and direction of the real photon permits t to 
be calculated without using the measured energy of the real photon, which is the quantity 
subject to the largest uncertainty. Thus, the value of t is calculated as 



COS I 



7 7^ 



1 + ifc - V^^' + COS0^*^) 



(3.1) 



for the exclusive event sample. The quantity —t is required to be smaller than 0.7 GeV^. 
The error caused by applying this expression to inelastic events is accounted for in the MC 
simulation that is used to calculate the fractional contribution of background processes per 
kinematic bin in xb, Q^, and —t. 

4. Extraction of asymmetry amplitudes 

The experimental yield M can be parameterized as 

AA(e,,P,,0) =/:(e,,P,)r?(ef,0)auu(0) x [l + P^A^^^^^c^) + eePiAi^icj)) + eeAc{(p)] ■ 

(4.1) 
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Here, C, is the integrated luminosity, r/ the detection efficiency, the beam polarization 
and (Tuu the cross section for an unpolarized target averaged over both beam charges and 
helicities. The asymmetries ^lu(^)' •^lu^^('!^)' ^'^'^ ^c(</') relate to the Fourier coefficients 
in Eqs. |1.2| - |1.4| according to Eqs. |2.2| - |2.4| . They are expanded in <f) as 



Am - E 4u5l'^ sin(-</') + E <u!r cos(n</)), (4.2) 

n=l n=0 

- E ^Swcs sin(n<^) + E <uSt^cs cos(n0), (4.3) 



ra=l n=0 



-4c(</') - E ^c"^""^^ cos(n0) sin (4.4) 

n=0 

where the approximation is due to the truncation of the in general infinite Fourier series 



caused by the azimuthal dependences in the denominators of Eqs. 2.2-2.4. The asymme- 



try amplitudes ^lu*!' ^c""'^ ^^°>'(o<^) relate to the twist-2 Fourier coefficients of the 



interference term appearing in Eq. |1.4| and further developed in Eqs. 1.5, 1.6 and 1 



respectively. The asymmetry amplitudes A^u'^vcS' ^Lu'^i''^''' ^c^^^"^^ ^"^^ 



lated to other Fourier coefficients in Eqs. |1.3| and 1.4, which are also explained above. The 
remaining asymmetry amplitudes are expected to be zero but were introduced to test for 
instrumentally induced harmonics. (The asymmetry amplitude ^lu DVCS ™ addition. 



arise through the interplay of numerator and denominator in Eq. 2.3 if the twist-3 Fourier 



coefficient sP^*-'^ has a sizeable value.) Comparison of predictions based on GPD models to 



data for either asymmetry amplitudes or the Fourier coefficients in Eqs. 1.3 and 1.4 provides 
similar information, as the asymmetry amplitudes relate to the corresponding Fourier co- 
efficients in an approximately model-independent way. This is due to the fact that the BH 
coefficients and lepton propagators are precisely calculable, and that in the kinematic 
region of HERMES the contributions from the squared DVCS term in the denominators 



of Eqs. 2.2-2.4 are expected to be much smaller than the BH contributions [pq ]. 

The extraction of the asymmetry amplitudes in each kinematic bin of xb, and t is 
based on the maximum likelihood technique |2^, providing a bin-free fit in <^ to the data. 
Its application here is explained in detail in Ref. [^5[. In the fit, weights were introduced to 
account for luminosity imbalances with respect to beam charge and polarization. No bal- 
ancing procedure was required for the target polarization, as the time averaged polarization 
was negligible as mentioned above. 

5. Background corrections and systematic uncertainties 

Each extracted asymmetry amplitude A is corrected, in each kinematic bin, for the contri- 
bution of semi-inclusive and exclusive background, which is mostly due to the production 
of vr*^ and r\ mesons. The background-corrected amplitude is calculated as 

. A /semi^semi /excl^excl /'c; i 

^corr — '-. 7 7 • I'^.ij 

-L /semi Jexcl 
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The fraction /semi of semi-inclusive background per bin is calculated using a MC simu- 
lation described in detail in Ref |^. It varies between 0.6% and 12.7% depending on 
the kinematic bin, with an average value of 3.3%. Based on the model in Ref. |3^, the 
fraction /exci of exclusive background is estimated in each kinematic bin and found to be 
below 0.7% everywhere. A direct search for exclusive neutral pions in the HERMES data 
supports this estimate . The semi-inclusive (exclusive) background can have a non-zero 
asymmetry Agemi (^exci)) as has been measured for, e.g., the semi-inclusive production 
of vr'^ mesons, which exhibits a sinusoidal (j) dependence on the beam helicity (3^. The 
beam-charge-dependent background asymmetry is zero at leading order QED. Hence the 
contributions from semi-inclusive and exclusive background constitute a dilution of Ac and 
effectively also of A\^u, as the latter involves only charge differences. In order to correct 
-^LU*^^ for the semi-inclusive background, the size of the beam-helicity asymmetry for this 
background is extracted from data by reconstructing neutral pions with a large fractional 
energy E^o/v > 0.8 (see Ref. [^] for details). For the exclusive background, the asymmetry 
cannot be extracted from data due to the small yield of exclusive pions. As the asymmetry 
is in the range [-1,1], a value of zero is assumed with a "statistical" uncertainty of 2/-v/T2, 
i.e., one standard deviation for a uniform distribution. The statistical uncertainty on the 
total correction due to the statistical uncertainties in the background fractions and asym- 
metries is propagated as a contribution to the final statistical uncertainty. In addition, 
half the size of the actual correction is assigned as systematic uncertainty to account for 
assumptions and approximations inherent in the approach. 

The dominant systematic uncertainty for the amplitude ^lu'dvcs ^° back- 

ground correction. In the case of ^"'^^ and ^ systematic uncertainty is pre- 

dominantly due to the combined contributions of possible deviations of the detector and/or 
the beam from their nominal positions ('alignment'), detector acceptance including smear- 
ing, and finite bin width in XB,t and Q^. The systematic uncertainty arising from these 
combined contributions is estimated by MC simulations using the GPD model described 



in Ref |33]. Note that a mistake has been found in this GPD model [^4[. However, the 



model of Ref. |33] described HERMES beam-charge |25] and preliminary (single charge 



beam-helicity asymmetries [35| and thus is considered adequate for systematic studies. In 
each bin, the systematic uncertainty is taken as the difference between the model predic- 
tions at the mean kinematic value of that bin and the respective amplitude extracted from 
the reconstructed MC data. 

Further systematic uncertainties arise from an observed relative shift of the squared 
missing-mass spectra between the electron and positron sample, with a magnitude of ap- 
proximately 0.2 GeV^ 1^^. The boundaries defining the exclusive sample in the missing 
mass spectra were adjusted to account for this shift. One quarter of the effect on the 
extracted asymmetries is applied as systematic uncertainty. The impact of both trigger 
and tracking inefficiencies was studied and found to be negligible. Also not included is any 
contribution due to additional QED vertices, as the most significant of these has been esti- 
mated to be negligible in the case of helicity asymmetries |37|. The systematic uncertainties 
are added in quadrature and given in Table ||. 
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Amplitude 


A lb (5(stat.) ± i5(syst.) 


Missing mass shift 


Background corr. 


Alignment, acceptance, bin width 


A sin 
^LU,I 


-0.224 ±0.028 ±0.020 


0.001 


0.005 


0.019 


^sin0 
^LU.DVCS 


0.043 ± 0.028 ± 0.004 


0.000 


0.004 


0.000 


.cos(O0) 


-0.020 ± 0.006 ± 0.008 


0.004 


0.000 


0.007 


* COS (j) 


0.055 ± 0.009 ± 0.004 


0.001 


0.001 


0.003 



Table 1: Overall values of the asymmetry amplitudes of greatest interest, at average kinematics 
{-t) = 0.12 GeV^, (xb) = 0.09, and (Q^) ^ 2.37 GeV^ The rightmost three columns present 
dominant contributions to the systematic uncertainties. Not included is a 2.8% scale uncertainty 
of the beam-helicity asymmetries due to the beam polarization measurement. 



6. Results 

Table |l| presents the asymmetry amplitudes of greatest interest extracted in the entire 
HERMES kinematic acceptance ("overall" results). The sm(j) amplitude of the beam- 
helicity asymmetry sensitive to the interference term is shown in the first row of Fig. |2|. 
It exhibits a large overall value of ^luI = -0.224 ± 0.028(stat.) ± 0.020(sys.), with no 
significant dependence on any of the kinematic variables —t, xb, and Q^. This implies a 
rather strong dependence of this amplitude on t for smaller values of —t, as the asymmetry 
amplitude has to vanish in the limit of vanishing —t due to the vanishing factor k in Eq. |1.5| . 
(In the limit of vanishing t, Cq^ remains finite and the dependences of K^h, Ki, Vi{(j)) and 
7^2(0) on t cancel in Eq. |2.2| .) The sin (j) amplitude of the beam-helicity asymmetry sensitive 
to the squared DVCS term is shown in the second row of Fig. ^. It also shows no kinematic 
dependence, with an overall value of ^LutDVCS — 0.043 ± 0.028(stat.) it 0.004(sys.). As 
explained above (see Eq. |2.1| ), the beam-helicity asymmetries measured previously with a 
single beam-charge are sensitive only to the combination of the results presented here, i.e., 
the single-beam-charge results are given as 

^LU l^^J ~ ^LU,I ^LU,DVCS' \^-^) 
if the contributions cjj from the interference term in the denominator of Eq. 2.1 can be 



neglected. Previous HERMES measurements [p4| , 25 1 found these contributions to be small 
compared to the remainder of the denominator, and a more precise constraint is presented 
below. Using the present data, the separate analysis of the positron {electron} data yields 
values for A'i'lf{ei) of -0.177 ± 0.022(stat.) {0.255 ± 0.051(stat.)}, in agreement with 
-0.181 ±0.046(stat.) {0.267 ± 0.065(stat.)} calculated from Eq. |6l| neglecting correlations 
from the commonality of the data sets. 

The sin 2(j) amplitude of ^lu,i is shown in the third row of Fig. ^. It has an overall value 
consistent with zero (—0.035 it 0.028 it 0.002) and thus, as expected, is much smaller than 



the corresponding sine/) amplitude. Those asymmetry amplitudes included in Eqs. 4.2-4.4 
as tests for instrumental effects are found to be compatible with zero. 

The data in Fig. |2| are compared with theoretical calculations to leading order in QED 
and QCD. In the GPD based model (VGG) described in Refs. |3^, 35], the dependences on 



^ and t are factorized while those on x and t may be entangled when a Regge-motivated 
ansatz is invoked. The model is formulated as a double distribution jl], ^ complemented by 
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Figure 2: The first (second) row shows the s'm(j> ampUtude of the beam-helicity asymmetry Altj.i 
(^LU,DVCs)i which is sensitive to the interference term (squared DVCS term), extracted from the 
1996-2005 hydrogen data in the entire experimental acceptance, and as a function of — i, a;B, and Q^. 
The third row shows the sm2(j) amplitude of Alu,i- The error bars (bands) represent the statistical 
(systematic) uncertainties. Not included is a 2.8% scale uncertainty due to the beam polarization 
measurement. The calculations arc based on the recently corrected minimal implementation |3^ 
of a dual-parameterization GPD model (Dual-GT) and on a GPD model ^0, 38 based on double- 
distributions (VGG). Both models use a Regge-motivated i-dependence. The band for the VGG 
model results from varying the parameters 5vai and bsca between unity and infinity. The bottom row 
shows the fractional contribution of associated BH production as obtained from a MC simulation. 



a D-term 1 39 ] , where the kernel of the double distribution contains a profile function 1 4C , 41 1 
that determines the dependence on ^, controlled by a parameter b ||4^. In the limit 6 — > oo 
the GPD is independent of ^. Note that 6^^! (^sca) is a free parameter for the valence 
(sea) quarks and thus can be used as a fit parameter in the extraction of GPDs from 
hard-electroproduction data. 

In each kinematic bin, a range of theoretical predictions was calculated [43| by varying 
the model parameters of only the GPD H, since these data are sensitive mostly to this 
GPD as explained above. Variants of the model are distinguished by differences in the t 
dependence of the GPD H, for which either a simple ansatz is used where the t dependence 
factorizes from the dependence on the other kinematic variables, or the Regge-motivated 
ansatz is employed. Since the differences are found to be small for all amplitudes shown 
in Fig. |2|, only the results based on the latter ansatz (VGG Regge) are displayed. The 
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broad width of the bands is due to the fact that the parameters 6vai and 6sea were varied 
between unity and infinity, with the variation in 6gea having the strongest effect. Note that 
including or neglecting the D-term in the GPD model does not change the result since 
it contributes only to the real part of the DVCS amplitude. The other model presented 
here (Dual-GT) is based on the corrected |34| minimal implementation [^] of the dual 
parameterization GPD model |44|, in which the dependence on ^ is factorized from the 
dependences on x and t. The t dependence in this model is also Regge-motivated. All 
models overestimate the magnitude of ^Lut approximately a factor of two. They are 
consistent with the observed shapes of the kinematic dependences onxB and (but not t). 
The part of the VGG band closest to the data, i.e., with the smallest absolute amplitude, 
corresponds to &sca = oo. Note that these models are for the elastic part of the cross 
section only while this measurement includes associated production in which the nucleon 
in the final state is excited to a resonant state. In the following it is considered whether 
the contribution from associated production can account for the observed discrepancies 
between model predictions and data. 

The bottom row of Fig. |2| shows the estimated fractional contribution from associ- 
ated BH production in each kinematic bin, calculated using MC simulations described in 
Ref. [^ 1. The overall value is about 12%. In an attempt to estimate ^Lut ™ elastic 
and associated production separately, the strong dependence of the fractional contribu- 
tions of elastic and associated production on the missing mass value in the exclusive region 
— (1.5 GeV)^ < < (1.7 GeV)^ can be utilized. The exclusive region can be split in sev- 
eral bins, each bin with its background-corrected amplitude ^corr = /eias^eias + /asso^asso- 
The fraction /das (/asso) of elastic (associated) production per bin is taken from the MC 
simulation, in which the DVCS process is not implemented because the BH cross section 
is expected to dominate that of DVCS not only for elastic but also for associated pro- 
duction Assuming that the values of A^u'j for elastic and associated production do 
not depend on M^, the two unknown amplitudes Aeias and j4asso are extracted from the 
five equations corresponding to the five M^-bins. The resulting overall sin0 amplitude 
from elastic production is found to be —0.209 it 0.066 and thus hardly differs from that 
reported in Table 0, while the one from associated production can only be constrained to 
be between —0.68 and 0.09 within one standard deviation in the statistical uncertainty. 
According to theoretical calculations a correction factor of 1.1 has to be applied to 
the measured beam-helicity asymmetry in HERMES kinematics due to the A resonance 
region {W < 1.35 GeV). For the associated DVCS amplitude these calculations are based 
on a model for transition GPDs, which are related to those on the nucleon within this 
model. Thus neither the extracted sin (f) amplitudes ^oias and ^asso nor model calculations 
support extreme scenarios in which the sincj) amplitude from associated production has 
an overall value of unity, which would be required to obtain a siiKp amplitude for elastic 
production as large as the value —0.39 or more predicted by the models shown. Thus 
associated production cannot account for the overestimate of the asymmetry amplitudes 
by the models. 

A promising alternative to comparing the data with existing models is to use a flexible 
GPD parameterization and perform a global fit to all DVCS data. First steps in this 
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direction have been made Pq , |47| , 48, 49, one of which found that a prehminary 
HERMES result on the sin (jj amphtude of the beam hehcity asymmetry for a single charge 
(positron) can be described by a fit to other DVCS data |Q. In order to provide 
additional input for future fits, in particular for the entangled ^ and —t dependences of 
GPDs, the amplitudes already presented in Fig. are shown in Fig. as a function of —t 
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Figure 3: The first (second) row shows the sin0 amplitude of the beam-helicity asymmetry Alu,i 
(^LU.DVCs) sensitive to the interference term (squared DVCS term), extracted from the 1996- 
2005 hydrogen data as a function of — i for three xb ranges. Correspondingly, the third row 
shows the sin(20) amplitude of Alu.i- The error bars (bands) represent the statistical (systematic) 
uncertainties. Not included is a 2.8% scale uncertainty due to the beam polarization measurement. 

for three different ranges of x-q. The possibly negative sva.2(j) amplitude for the largest x-q 
bins in Fig. |^ is found to be independent of t in the lower right panel of Fig. |3|. 

The cos{n(j)) amplitudes (n = 0-3) of the beam charge asymmetry are shown in Fig. |^. 
The cos(0(/)) and the coscf) amplitudes, i.e., the amplitudes related to twist-2 GPDs, are 
zero at small values of —t and become non-zero with increasing values of — t, with opposite 
sign and smaller magnitude for cos(0(?l)) as expected from Eq. |1.8| . It is interesting to note 
that A^Q^'^ and ^lu'i show a fundamentally different dependence on — t, despite relating 
to the real and imaginary parts of the twist-2 helicity-conserving DVCS amplitude via the 
same factor k oc \f—i/Q in Eq. and Eq. |L^, respectively. The coscj) amplitude does 
not exhibit any kinematic dependence on x-q or Q^. It is in agreement with the earlier 
HERMES measurements based on subsamples of the data used in this analysis p5|| . 
The cos(20) amplitude, which is related to twist-3 GPDs, is suppressed as expected and 
found to be compatible with zero. Also, the cos(30) amplitude, which is related to gluon 
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Figure 4: The cos(n(/)) amplitude {n = 0-3) of tlie beam-charge asymmetry Ac, extracted from 
the 1996-2005 hydrogen data in the entire experimental acceptance, and as a function of —t, xb, 
and Q^. The error bars (bands) represent the statistical (systematic) uncertainties. The theoretical 
calculations are based on the models that are unable to describe the data in Fig. ||. For the VGG 
model the parameter settings 6vai = oo and fegca = 1 are used and the contribution from the D-term 
is set to zero. The bottom row shows the fractional contribution of associated BH production as 
obtained from a MC simulation. 



helicity-flip GPDs, is found to be consistent with zero. No striking additional features are 
observed in Fig. ^ where the cos{n(f)) amplitudes are shown as a function of —t for three 
distinct xb ranges. 

The theoretical calculations shown in Fig. |^ are based on either the Dual-GT or the 
VGG model. For the VGG model the parameter settings 6vai = oo and 6sca = 1 are used 
and the contribution from the D-term is set to zero, as only this set of parameters yields 
a good description of the BCA data [p^ , |25|| . Note that the same set, in particular the 
setting bsea = Ij leads to amplitudes with the largest magnitude among those represented 
in the bands in the top row of Fig. ^, i.e., it clearly does not describe the data related to 
the imaginary part of the DVCS amplitude. It appears that additional degrees of freedom 
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Figure 5: The cos(n0) amplitude (rt = 0-3) of the beam-charge asymmetry Ac, extracted from the 
1996-2005 hydrogen data as a function of —t for three a;B ranges. The error bars (bands) represent 
the statistical (systematic) uncertainties. 

in the calculation of the BCA, such as the value assigned to the D-term, allow the VGG 
model to be tuned to resemble the BCA data. Similarly, the Dual-GT model does not 
describe the data in Fig. ^ but is in reasonable agreement with the BCA data shown in 
Fig. ^. (The sudden increase of the cos (f> amplitude predicted by this model in the highest 
xb and bins might be due to the fact that this model is designed for small and medium 
values of xb up to 0.2.) While the increase {decrease} of the cos cp {cos(0(/))} amplitude 
with —t is well reproduced within these models, the contribution of associated processes 
not included in these models is expected to also increase with —t as shown in the bottom 
row. 

7. Summary 

Previously unmeasured charge-difference and charge-averaged beam-helicity asymmetries 
in hard electroproduction of real photons from an unpolarized proton target are extracted 
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from data taken with electron and positron beams. The sm(j) amphtudes of these beam- 
hehcity asymmetries are sensitive to the interference term (twist-2) and to the squared 
DVCS term (twist-3), respectively, whereas earlier measured beam-helicity asymmetries 
with a single beam-charge are sensitive to only their linear combination. In addition, the 
most precise determination of the beam-charge asymmetry is presented, which provides 
access to the real part of the DVCS amplitude. The GPD models presented are not able to 
describe the sin amplitude sensitive to the interference term, while they can be adjusted 
to resemble the results on the beam-charge asymmetry, presumably because the model 
calculations have additional degrees of freedom in the latter case. The amplitudes related 
to higher-twist or gluon helicity-fiip GPDs are found to be compatible with zero. The 
results presented on these charge-decomposed beam-helicity asymmetries and on the high- 
precision beam-charge asymmetry have the potential to considerably constrain the GPD 
H when used in comparison with future GPD models or as input to global fits. 
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A. Tables of Results 
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Table 2: Bin sizes, average kinematic values and results of the asymmetry amplitudes presented 
in Fig. I 
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Table 3: Bin sizes, average kinematic values and results of the asymmetry amplitudes presented 
in Fig. |. 
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Table 4: Bin sizes, average kinematic values and results of the asymmetry amplitudes presented 
in Fig. |. 
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Table 5: Bin sizes, average kinematic values and results of the asymmetry amplitudes presented 
in Fig. ||. 
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